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ABSTRACT 

We have found a system listed in the Kepler Binary Catalog (P or b = 3.273 days; Prsa et al. 2010) that we 
have determined is comprised of a low-mass, thermally-bloated, hot white dwarf orbiting an A star of about 
2.3 M Q . In this work we designate the object, KIC 10657664, simply as "KHWD3" ("Kepler Hot White Dwarf 
3"). We use the transit depth of ^0.66%, the eclipse depth of ~1.9%, and regular smooth periodic variations 
at the orbital frequency and twice the orbital frequency to analyze the system parameters. The smooth periodic 
variations are identified with the classical ellipsoidal light variation ("ELV") and illumination ("ILL") effects, 
and the newly utilized Doppler boosting ("DB") effect. Given the measured values of R/a and inclination angle 
of the binary, both the ELV and DB effects are mostly sensitive to the mass ratio, q = Mz/M\, of the binary. 
The two effects yield values of q which are somewhat inconsistent - presumably due to unidentified systematic 
effects - but which nonetheless provide a quite useful set of possibilities for the mass of the white dwarf (either 
0.26±0.04M Q or 0.37 ± 0.08 Mq). All of the other system parameters are determined fairly robustly. In 
particular, we show that the white dwarf has a radius of 0.15 ±0.01/?© which is extremely bloated over the 
radius it would have as a fully degenerate object, and an effective temperature T e ff ~ 14, 500 K. Binary evolution 
scenarios and models for this system are discussed. We suggest that the progenitor binary was comprised of a 
primary of mass ^2.2 M Q (the progenitor of the current hot white dwarf) and a secondary of mass ^1.4 M Q 
(the progenitor of the current A star in the system). We compare this new system with three other white dwarfs 
in binaries that likely were formed via stable Roche-lobe overflow (KOI-74, KOI-81, and the inner Regulus 
binary). 

Subject headings: techniques: photometric — stars: binaries: eclipsing — stars: binaries: general — stars: 
evolution — stars: variables: other — stars: subdwarfs 



1. INTRODUCTION 

The exquisite photometric precision of the Kepler mission 
has led, among many other things, to the discovery of two 
systems in which a hot white dwarf both transits its parent 
star, and is eclipsed by it (Rowe et al. 2010; van Kerkwijk 
et al. 2010). In the case of KOI-74, the relative transit and 
eclipse depths are only 5 x 10 -4 and 12 x 10~ 4 , respectively, 
while the amplitudes of the ellipsoidal light variations and 
Doppler boosting effect have even more remarkably small am- 
plitudes of 1 .4 x 10 -4 and 1 .0 x 10 -4 , respectively. Such small 
effects are unexplored with ground-based astronomy. 

The discovery of transiting/eclipsing white dwarfs which 
exhibit the effects mentioned above are amenable to many in- 
teresting system diagnostics, and can thereby be potentially 
very revealing about the binary stellar evolution that has led 
to their current configuration (e.g., van Kerkwijk et al. 2010; 
Di Stefano 2010). In particular, if the progenitor of the white 
dwarf was initially of mass < 2.2M© the orbital period may 
be tightly correlated with the white dwarf mass (e.g., Rap- 
paport et al. 1995). In all cases where the current-epoch bi- 
nary orbital period is short (e.g., less than a month) the sys- 
tem almost certainly involved a phase of mass transfer from 
the progenitor of the white dwarf to what is presently the nor- 
mal stellar companion. In the process, a substantial amount of 
the original primary mass was transferred to its companion or 
lost from the system. The temperature of the white dwarf pro- 
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vides some important information on the cooling history and 
the age since the mass transfer event (see, e.g., van Kerkwijk 
et al. 2010; Di Stefano 2010). 

In this work, we report the discovery of a third hot white 
dwarf (with T e ff ~ 14,500 K) orbiting an A star, which we 
herein designate as "KHWD3". The star, KIC 10657664, was 
identified and catalogued by Prsa et al. (2010) in the Kepler 
binary star catalog, although it was not identified as harboring 
a hot white dwarf companion. Its transit and eclipse depths 
are about an order of magnitude larger than those of KOI-74 
and KOI-8 1 , due mostly to the fact that the radius of the white 
dwarf in KHWD3 is a factor of >3 times larger than in the 
systems discovered earlier. 

The layout of this paper is as follows. In §2, we present 
the Kepler light curve for this system, and the analyses which 
yield the flux ratio of the two stars, the ratio of their radii, the 
orbital inclination angle, and the first three harmonics of the 
light curve which reveal ellipsoidal light variations, an illumi- 
nation effect, and Doppler boosting. In §3 we determine the 
model-independent parameters of the binary system includ- 
ing the mass, effective temperature, radius, and luminosity 
of the white dwarf, and the mass ratio. Inferences about the 
properties of the primary star allow us to estimate its mass 
at 2.3±0.4M Q . From this estimate and the measured mass 
ratio we arrive in §4 at two possible solutions for the mass 
of the white dwarf: either 0.26 ± 0.04 M Q or 0.37 ± 0.08 M Q , 
depending on whether we emphasize in our analysis the ellip- 
soidal light variations or the Doppler boosting effect. In §5 
we discuss various evolutionary scenarios that can lead to the 
present system. 

2. LIGHT CURVE ANALYSIS 

The uncorrected light curve for KHWD3, downloaded from 
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the Kepler public archival is shown in Figure [T] The data 
are comprised of ~ 10 days of quarter zero (Q0) data and 
~ 33 days of quarter one (Ql) data^. Outside of the easily- 
identified eclipses, the variation of the light curve is domi- 
nated by a low frequency trend, assumed to be a systematic 
associated with the Kepler spacecraft. The remaining periodic 
variation on times at or shorter than the orbital period and in 
phase with the conjunctions is assumed to be astrophysical 
in nature (as a result of, for example, ellipsoidal light varia- 
tions and/or Doppler boosting). We accounted for these three 
effects in a light curve model, as described below, and deter- 
mined parameters of the eclipsing system and the harmonic 
content of the out-of-eclipse variation. 

2.1. Light Curve Model 

We modeled the eclipse events assuming two spherical stars 
having radius ratio R2/R1, and observed flux ratio FijF\, as- 
signing the indices of 1 and 2 to primary and secondary, re- 
spectively. The binary was constrained to a Keplerian orbit 
parameterized by a period P, a normalized semi-major axis 
distance a/R\, an inclination to the sky plane i, an eccen- 
tricity e, an argument of periastron ES and a time of inferior 
conjunction (i.e., at the time of the eclipse of the primary), f; c . 

The normalized eclipse light curve, /(f), was calculated to 

be 
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where z(t) is the sky-projected separation of the centers of the 
two stellar components and A is the fraction of the eclipsed 
disk blocked by the occulter, given analytically by Mandel 
& Agol (2002). The limb darkening coefficients u and v pa- 
rameterize the radial brightness profile, I(r), of either binary 
component as 
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where r is the projected radial distance from the stellar center, 
normalized to 

The out-of-eclipse light curve was modeled as the product 
of a harmonic series, h(t), and a low order, quarter specific 
polynomial, pnit) (with Q = Q, 1); 



h(t) = 1 + ^A k sin [k<j)(t)] + Y, B k cos [*#»)] (3) 

k=l k=l 

Nq 

PQ(t)=J^c®Ht-t y (4) 

for a set of amplitudes (A^, B^) and polynomial coefficients 
[CJP] where 4>{t) = (27r/.P)(f-fj c ), and fo is some fixed time 
near the start of Q0, No = 3 and Ni=6. 

The full light curve model, 1(f), is given as the product of 
all three variations, i.e., £(t) = f(t)h(t)p Q (t). 

We fitted the full model to the data with R2/R1, a/R\, 
F2/F1, wi, U2, P, i, e, 03, fie, Ak, Bk and C, ( , e ' as free param- 

5 \protect ittp://archive. stsci.edu/kepler/publiclightcurves. html 

6 Kepler "quarters" are continuous observing blocks separated by short 
times in which data are downloaded from the spacecraft. Only Q0 and Ql 
are currently available to the public. We have utilized only the unprocessed 
("raw") data products in our analysis. 



eters. We fixed the quadratic limb darkening coefficient of 
the primary to vi = 0.2964, as would be expected (according 
to Sing 2010) for a star observed in the Kepler bandpass and 
having the stellar parameters T e ff = 10000 K, \ogg = 4.0 and 
[M/H] = -0.10 as determined for KHWD3 and tabulated in 
the Kepler Input Catalog. We adopted a linear limb darken- 
ing model for the secondary by fixing the quadratic term V2 to 
zero; the inclusion of this degree of freedom had little effect 
on the results but was added for completeness. 

The continuously defined model, £(t), was numerically in- 
tegrated before being compared with the long cadence Ke- 
pler light curve. In detail, for each measured time, tj, we 
take iij uniform samples tj^ = tj+kAtj — T^t/2, separated by 
Alj = T int/"> over the long cadence integration interval of 
ri n t = 30 minutes. Then, these times were evaluated with l(t ) 
and averaged to determine the integrated model flux, Tj, at 
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The required number of samples, rij, depends on the curvature 
of the model relative to the photometric precision at a given 
time; iij > 1 was only necessary during the relatively sharp 
eclipse and transit events. In practice, we determined tij by 
increasing its value until the goodness-of-fit statistic (defined 
below) plateaued. We settled upon rij = 60 for any tj within 
~ 2 hr of an eclipse or transit and rij = 1 elsewhere. 

We determined the best fit model to the data by minimizing 
the x 2 goodness-of-fit statistic, defined as 



(T s -F s ) 2 



(6) 



where F s is the measured flux at time t s and a is the expected 
statistical error in the flux measurements. We selected a = 
130 ppm such that the reduced x 2 was unity for the best fit 
solution. 

The solid curves in Figures Q] [2] and [3] trace the best fit 
solutions to the data. 

2.2. Markov chain Monte Carlo Analysis 

We determined the posterior probability distributions for 
the fitted parameters by employing a Markov chain Monte 
Carlo (MCMC) algorithm. With the MCMC algorithm, first, a 
randomly chosen parameter is perturbed by an amount drawn 
from a normal distribution of some fixed width. Next, this 
new parameter and the remaining unperturbed parameters are 
used to compute x 2 t as defined in Eqn. (0] to determine 
a likelihood £ cx exp(~x 2 /2). The new parameter is "ac- 
cepted" according to a Metropolis-Hasting jump condition: If 
the likelihood £ is greater than the unperturbed likelihood C 
then the perturbation is accepted, otherwise it is accepted with 
a probability of £ / C. After many perturbations, the result- 
ing "chain" of accepted parameters samples the desired joint 
posterior probability distribution. For a more detailed intro- 
duction to the MCMC algorithm, refer to the appendix of the 
work by Tegmark et al. (2004). 

We generated a chain of 2 x 10 6 links, having selected per- 
turbation widths such that ~ 40% of jumps were accepted. 
The chain was checked for adequate mixing and convergence 
by visual inspection, and by observing that the number of 
links was much larger than the autocorrelation length (equal to 
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FIG. 1 . — Kepler light curve for KIC 10657664. The solid line gives the best fit solution to the data for the model given in §2. The dashed curve is the estimated 
systematic component of the full model. The lower panel shows residuals from subtracting the best fit model from the data (in units of parts per million). 
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FIG. 2. — Period-folded, systematic corrected and normalized Kepler light curve for KIC 10657664 near the phases of eclipse and transit. The solid lines give 
the best fit solution for the model described in §2. The lower panel shows residuals from subtracting the best fit model from the data. 



the number of links at which the chain autocorrelation drops 
below one half) for any selected parameter. We report, in Ta- 
ble [JJ the 50%, 15.8% and 84.2% values (corresponding to 
the median and ±l-sigma values) of the cumulative distri- 
bution for each parameter, marginalizing over the remaining 
parameters. In the upper-left corner of Figure|4j we show the 
joint distributions for the parameters R2/R1, Fz/F\, a/R\, i 
and Mi to demonstrate their correlations. The remaining pa- 
rameters, including the harmonic function amplitudes, were 
weakly correlated with other parameters. 

2.2.1. Out-of-eclipse harmonic content 

We measured amplitudes (being significantly different from 
zero) for, in order of decreasing amplitude, the cos(20), 
cos(</>), sin(</>), sin(2^>) and cos(30) modes of the out-of- 
eclipse variation. Repeating the above analysis including only 
an integral number of orbital cycles did not modify these re- 
sults. The four largest modes are plotted individually for com- 
parison in Figure [3] The fitted values of the 6 harmonic am- 
plitudes are summarized in Table 1 . 



In §3, we associate these modes with three possible phys- 
ical effects and describe an analysis used to determine addi- 
tional system parameters (such as the mass ratio q), assuming 
the validity of these associations. 

Correlated Noise - We have assumed a "white" noise spec- 
trum in both our best-fit and MCMC analyses. However, we 
note that correlated noise could have a significant affect on 
our estimates of harmonic content. To investigate this possi- 
bility further, we completed a secondary analysis, as follows, 
whereby we determined what we consider to be the most con- 
servative errors on these amplitudes. 

We first corrected for the low-frequency trend by divid- 
ing the data by the best-fit estimate for po(t). Then, having 
fixed the best-fit estimates for the eclipse parameters, we fit a 
unique harmonic model to each individual cycle (i.e., by using 
only the data from one primary transit to the next). We then 
determined the mean amplitude, for each mode, and the stan- 
dard deviation in that mean averaging over the 12 available 
cycles. 

We found that the scatter across all cycles, in a given har- 
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monic amplitude, was consistent with random noise (i.e., no 
minority of cycles was significantly biasing the mean). The 
means across all modes agree with the MCMC-determined 
amplitude medians within the MCMC-determined errors. The 
standard deviations in the means were larger, as may be ex- 
pected in the presence of correlated noise, than the amplitude 
uncertainties determined via MCMC. However, only the sin <j> 
and cos <f> modes have substantially higher uncertainties, with 
each being roughly twice as large as the MCMC estimates. 
These inflated, and likely overestimated, errors translate into 
a ~ 10% correction in the uncertainty in the secondary mass 
estimate (as determined in §3). All other system parameters 
were less significantly affected. We therefore decided to use 
the MCMC error estimates in the subsequent analysis. 

2.2.2. Eccentricity 

We also report a weak detection of eccentricity for this 
system, having measured a non-zero value for ecos(03) = 
0.0029 ± 0.0005 as inferred from the — 4.5 minute delay (rel- 
ative to the circular orbit expectation) of the start of secondary 
eclipse. We note that we did not account for delays as a re- 
sult of the finite speed of light crossing the orbit; however, 
given the estimated separation of the binary (see §3), we ex- 
pect light travel time delays of ~ 30 seconds. This correction 
corresponds to an error in ecostt! that is smaller than the esti- 
mated statistical error. In §3.6.1, we comment on the conse- 
quences of a non-zero eccentricity in the interpretation of the 
harmonic modes of the out-of-eclipse light curve. 

3. LIGHT CURVE MODEL-INDEPENDENT SYSTEM PARAMETERS 

From the measured values of R2/R1 and F>/F[ we may de- 
termine the mean surface brightness ratio Bz/B\ between the 
secondary and primary according to 

^ expjhc/XkTj)-! 
exp(/ic / \IcT2) - 1 

where in the second line we have substituted for the sur- 
face brightness ratio the expectation assuming black-body 
emission and a narrow observed wavelength range (with A « 



TABLE 1 

KHWD3 Light Curve Model Parameters 



Parameter 


Median 


84.2% (+lcr) 


15.8% (-lo-) 


P [day] 


3.273713 


+0.000008 


-0.000008 


fic [BJD] 


2454951.85857 


+0.00012 


-0.00012 


ecos(B5) 


0.0029 


+0.0005 


-0.0005 


|esin(G3)| 


$0.1 a 






Eclipse Parameters 








R 2 /Ri [xlO 2 ] 


8.100 


+0.024 


-0.020 


F 2 /F 1 [X10 2 ] 


1.900 


+0.002 


-0.002 


a/Ri 


7.02 


+0.08 


-0.10 


i [deg] 


84.46 


+0.14 


-0.18 


U\ 


0.20 


+0.02 


-0.02 


112 


unconstrained b 






Out-of-eclipse Harmonic Amplitudes [x 10 4 ] 




Ai[sin(0)] 


1.03 


+0.04 


-0.04 


B|[cos(<^)] 


2.64 


+0.04 


-0.04 


A 2 [sin(20)] 


-0.44 


+0.04 


-0.04 


B 2 [cos(20)] 


-5.68 


+0.04 


-0.04 


A 3 [sin(3<ffl 


0.02 


+0.04 


-0.04 


B 3 [cos(3<ffl 


-0.13 


+0.05 


-0.05 



" e sin 05 is nearly unconstrained by the data (and uncorrelated with the re- 
maining parameters). Here we report the full range of its weakly converged 
Markov chain. 

b The limb darkening parameter of the secondary was limited to positive val- 
ues less than 0.8. 



6000A for Kepler). To determine this brightness ratio in prac- 
tice, for a given set of temperatures, we integrated each black- 
body spectrum over the wide Kepler response function. 

The harmonic content of the out-of-eclipse light curve vari- 
ation for KHWD3 is presumed to be astrophysical, i.e., related 
to the binary system, rather than systematics. In the previous 
section, we reported on the measured values of the harmonic 
amplitudes, Ai and for k = 1, 2, or 3, where we remind the 
reader that the out of eclipse variation was modeled as 

3 3 
h(t ) = 1 + ]T A k sin [kcj)(t )] + B k cos N>(0] ■ (8) 
*=1 k=l 

Briefly, we expect that the sin(^) term is due to Doppler 
boosting (hereafter "DB"), the cos(0) term is largely due to 
mutual illumination (hereafter "ILL") between the two com- 
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ponents, and the cos(20) and cos(3</>) terms are largely due to 
the ellipsoidal shape deformation of the primary by the sec- 
ondary's tidal field (so-called ellipsoidal light variation, here- 
after "ELV"). In this nominal model, we do not expect non- 
zero amplitudes for the sin(2</>) or sin(30) modes. 

3.1. Expected harmonic amplitudes 

We assume a circular binary orbit throughout the following 
discussion. A departure from this assumption is discussed in 
§3.6.1. 

The predicted amplitudes At, B k may be decomposed into a 
sum of effect-specific amplitudes such that 



A k =A™+Af LV +A 1LL 



B k =Bf B +Bf LV + B 



k 

ILL 



(9) 
(10) 



In turn, the effect-specific amplitudes are a sum of the indi- 
vidual contributions from each member of the binary; for ex- 
ample,A° B =A A DB ' 1 +A° B ' 2 . 

Doppler boosting - The apparent spectral intensity of either 
component changes throughout the orbit as a result of Doppler 
shifting, photon emission rate modulation and beaming. The 
combination of these effects is called Doppler boosting; it is 
discussed in the context of Kepler by Loeb & Gaudi (2003) 
and by Zucker et al. (2007) and was detected for at least one 
of the two other candidate hot white dwarf companions in the 
Kepler field (KOI-74 and KOI-81; van Kerkwijk et al. 2010) 
and for one planetary system (CoRoT-3; Mazeh & Faigler, 
2010). The time variability of the Doppler boosting signal 
is carried by the projected radial velocity v,-n 2) where, for a 
circular orbit, v r (i2) oc sin</> and the total contribution to the 
harmonic out-of-eclipse content by Doppler boosting is 



(*)-(*) 



2na 
~Pc~ 



Oil 1 + 
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-a 2 (l+q) 



-1^2 

Fi 



(11) 



The DB prefactor (a) depends on the shape of the emission 
spectrum in the emitter's rest frame such that 



/C/ln/v 
v d\nv 



(12) 



where the average is weighted over the Kepler response. It 
has been empirically validated in the case of KOI-74 (Ehren- 
reich et al. 2010) that prefactors estimated assuming a black- 
body emission model are overestimated by greater than 20% 
for stars with T w 10000 K. As such, we utilize the synthetic 
spectra by Kurucz (1979) to generate DB prefactors over a 
range of effective temperatures, by fitting power law models 
(an impressively good approximation over the Kepler band 
pass). The fits are performed in a weighted least-squares sense 
with the weight at a given frequency being specified by the 
Kepler response function. For the A star primary, we selected 
model spectra of fixed gravity (log g = 4.0 [cgs]) over a range 
of temperatures from 8250 K to 1 1500 K and derived DB pref- 
actors. When then fitted a linear function in temperature to 
these prefactors for use in our subsequent analyses: 



ai(7i)« 1.90-1.12 



10000K 



(13) 



The DB prefactors for the WD component [a 2 (T 2 )] were 
computed analogously, but assuming a simple black-body 
emission model. 

Ellipsoidal light variation - The mutual gravitational inter- 
action between members of a binary will induce nominally 
prolate distortions of their surfaces having major axes along 
the line connecting their centers. As a result, the sky -projected 
cross-sectional area, and consequently the observed flux, are 
time variable (Morris 1985; see also Pfahl, Arras, & Paxton 
2008). For a circular orbit, the expected maxima of this vari- 
ation occur twice per orbit at the quadratures [generating a 
large cos(20) amplitude]. Kopal (1959) calculated the ELV 
contribution to the out-of-eclipse harmonic content as (where 
we follow the notation by Morris 1985): 



B 



ELV , 



-3Zi(3)<? 



jELV , 



5 



-Z\ (2)q ( — ) sin 2 
V a 

-5Z x (2>)q(— ) 



(14) 



(15) 



(16) 



where the prefactors Zj(2) and Zi(3) (Eqn. 3 in the work by 
Morris 1985) depend on limb-darkening and gravity darken- 
ing parameters with Zi(2) w 1.5 and Zi(3)/Zi(2) w 0.05. In 
our calculations, we assume a fixed linear limb-darkening pa- 
rameter [u = 0.44; according to Sing (2009) for the Kepler 
Input Catalog stellar parameters listed in §2.1] and a gravity 
darkening parameter that depends on temperature according 
to the von Zeipel law (Eqn. 10 in the work by Morris 1985). 

The contribution to the ELV due to the distortion of the 
secondary by the primary is smaller by a factor of at least 
(R 2 /R\f{F 2 /Fi)/q 2 ~ 0.05% and is therefore safely ignored. 

Mutual illumination - It is expected that radiation arising 
from one binary component and incident on the other will 
subsequently either be scattered or absorbed and reemitted. 
As such, the surfaces of the components that are directly fac- 
ing one another will be brighter than other regions. As the 
binary orbits, we perceive these phases as a modulation of 
the out-of-eclipse total flux level. The maximum illumina- 
tion of the primary is at inferior conjunction while that of the 
secondary is at superior conjunction, contributing with oppo- 
site signs to a predominately cos (j> mode of the out-of-eclipse 
variation. If the binary is in radiative equilibrium, then all in- 
cident radiation (bolometric) must be absorbed and reemitted 
at approximately the effective temperature of the illuminated 
member (Eddington 1926). In this scenario, Kopal (1959) 
again provides the expected out-of-eclipse harmonic content 
due to mutual illumination as 



1 16 V a 



17 (R 2 

16 \~~a~ 

17 (Ri 

16 V a 

17 /* 
16 \~a 



yILL , 



11/* 

3 4 V a 
11/* 
3 4 V a 

16 _3_ 

Hit 2 + 16 

16 3 
27tt 2 + 16 



TiVfRi^ 2 



BC 2 
BCi 

* 

a 

R 2 
a 



(17) 



h 
7] 
BC 2 
BC X 



(18) 



where BC is an approximate bolometric correction to the vi- 
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sual (see Allen 1964) with 



29,000K 

-2.51ogfiC = -42.5+101ogr eff + — . (19) 

?eff 

3.2. First cut at determining the constituent masses: an 
inconsistency 

In principle, we may solve for the system parameters M\ , 
7*i, R\,Mi, 7*), Ro and a using the formalism introduced in the 
previous sections (Eqns. 7, 11-19), Kepler's third law, and the 
results from the light curve analysis presented in §2. However, 
the illumination effect (given by Eqns. 17 and 18) is nearly in- 
dependent of the system masses, which are largely determined 
by the Doppler boosting (Eq. 1 1) and ELV effects (Eqns. 14- 
16). If we take the simplest forms of both the DB and ELV 
effects we find: 



:2.0 



/ 27Tfl\ 



V Pc J \l+q 

3 



q 



-1.5,(^1 

a 



(20) 
(21) 



where the coefficients a t and Zi(2) have been set equal to 
their nominal values, i has been set to 90° without loss of ac- 
curacy, and q = M%lM\. For purpose of this particular ex- 
ercise, we have neglected the relatively small effect of the 
Doppler boosting of the hot white dwarf. If we eliminate a 
in favor of P through Kepler's third law, and plug in our mea- 
sured values for P and Ri/a,we find: 



;9.6x 10" 4 to 



■4 1/3 



q 



(1+4) 2 / 3 



1.03 x 10" 



fif LV ~-4.3 x 10~ 3 q ~ -5.68 x 10" 4 



(22) 
(23) 



where the numerical values on the right sides of the equa- 
tions are taken from Table 1, and the mass of the primary, 
mi has been expressed in solar units. Since these equations 
depend on q and only very weakly on the actual mass of the 
primary, they both, in effect, independently determine q. For 
any primary mass between 2 and 3 M®, the two equations 
yield: q ~ 0.078-0.090 and q ~ 0.13, respectively, depend- 
ing on whether we utilize the DB or ELV amplitudes to de- 
termine q. These are clearly inconsistent, and this remains 
the case, even when we later utilize all the terms in Eqns. (7), 
(1 1)-(19), and take full account of the statistical uncertainties 
associated with the measure quantities. Since there is a de- 
pendence on Mi in one of the two equations, there exists a 
formal solution to the two equations for M\ and Mo (0.7 M Q 
and 0. 1 M®, respectively). However, this solution is not phys- 
ically plausible, and therefore we are left with this intriguing 
inconsistency in the value of q that we determine. 

3.3. Roughly determined parameters: Ti and pi 

Once we decide on the mass ratio for the system, we will 
need an estimate of the mass of the primary in order to deter- 
mine the mass of the white dwarf. There are two accessible 
parameters that we can utilize to estimate the primary mass: 
the effective temperature, 7] , and the mean stellar density, pi . 
T] is listed for the primary star in the Kepler Input Catalog 
(KIC) as 10,500 K, based on 5-color photometry. We have no 
readily available estimate of the uncertainty in the KIC value 
for 7i, but we estimate that it could be ~ ±500 K. As we shall 
see in §3.5, the illumination effect can be used to directly in- 
fer the effective temperature of the primary, and yields a result 



of 9500 ± 150 K, which would then be some 2 a away from 
the KIC tabulated value. We therefore consider a range for 71 
that includes both the KIC value and the one we derive from 
the illumination effect. 

The second readily, and more accurately, determined pa- 
rameter associated with the primary is its mean stellar density, 
pi = Mi/(4/3ttR]). We point out that pi is only a function 
of the light curve-determined parameters P and a/R\, with a 
weak dependence on q. This dependence of pi on these pa- 
rameters may be found by dividing Kepler's third law by R\: 



— ) =— GP 2 pi(\+q) 
RiJ 3tt H H 



(24) 



This relation was pointed out by Seager & Mallen-Ornelas 
(2003) in the context of transiting exoplanets. With exoplan- 
ets, the value of q is much smaller than the error in a/R\ such 
that pi can be determined independently of any mass infor- 
mation. While the same is not quite true with our binary, the 
value of q is still small enough (with q < 0.2 as a very conser- 
vative estimate) such that a moderately accurate value of pi 
may be determined independently of q. In particular, we find 
that pi fn 0.61 ±0.02 g cm by using Eqn. (24) with q = 
and the results from our photometric analysis. We consider 
this estimate to be maximal and accurate to within 20% of its 
true value. 

3.4. A prior on M\ 

In the previous section, we argued that 7j ~ 10,000 K and 
that pi < 0.61 g cm" 3 . Given these two values, we can moti- 
vate a prior on Mi . 

If we assume the primary is on the zero age main sequence 
(ZAMS), then we may utilize the approximate formula by 
Eggleton (2006) to estimate its mass as a function of either 
effective temperature or density. Picking the former to be 
in the range 9000 K < Ti < 11,000 K, we find a mass 2.0 
< Mi /M Q < 2.6. Choosing the latter to be in the range 0.5 1 g 
cm" 3 < pi < 0.61 g cm" 3 , we find a mass 2.5< Mi/Mq < 3.2. 
These ranges in Mi intersect, suggesting Mi ~ 2.5 M Q ; alter- 
natively, either the primary has aged somewhat off the ZAMS 
or the metallicity of the primary is different from Solar. To 
investigate these possibilities further, we utilized the Yonsei- 
Yale stellar evolution isochrones (Yi et al. 2001) to deter- 
mine a likely mass range subject to the priors on tempera- 
ture and density as given above and additionally a normally 
distributed prior on the metallicity, [Fe/H] = 0.0±0.2 (apply- 
ing the methodology described by Carter et al. 2009). From 
this analysis, we find that Mi = 2.3+Q4 M q with an estimated 
age for the star in the range 200-600 Myr; although, given a 
likely mass transfer history (see §5), the interpretation of this 
age estimate is unclear. 

Considering the above ranges in Mi , we opted for a conser- 
vative mass prior for the primary being normally distributed 
about Mi = 2.5 M with an rms width 0.5 M Q . 

3.5. Derived system parameters subject to the prior on Mi 

In the following analysis of the system parameters we 
assume that the illumination model is reasonably correct 
(Eqns. 17 and 18), and we utilize either the ELV effect 
(Eqns. 14-16) or the Doppler boosting effect (Eqns. 11 and 
12). Given the common prior on Mi (see the previous sec- 
tion), either effect (ELV or DB) results in an independent esti- 
mate of q. As discussed earlier, however, we find that these es- 
timates disagree with one another: ignoring DB [by nullifying 
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the statistical weight of the measurement of the sin(0) mode] 
gives i7~0.15i0.01, whereas ignoring ELV [by nullifying 
the statistical weight of the measurement of the cos(20) mode] 
gives g~0.11±0.01. Given this discrepancy, we opted to de- 
termine and report system parameters for these two scenarios 
separately. 

To determine posterior distributions for all system parame- 
ters, we executed another Markov chain Monte Carlo algo- 
rithm (as described in §2.2) subject to the likelihood £ oc 
exp(— x|/2) where 



(A|) 2 (a|) 2 (a-) 2 (A0 , (AP)2 



Rl/Rl 



F2/F1 



PR 2 /R\,a/R 



(A|)(A£) (a|) ( ao 



Pa/Ri 



OR^JR^a/R 



+ PR 2 /Ri 



&R 2 /R,Vi 



(AAfc) 



(MX-2.5M, 



Va/RxPi 
2 



+2. 

k 



a 



,-,) 



(O.5M ) 

where for any parameter "x", Ax is the difference between the 
median estimate of the parameter x and its current value (at a 
given link in the Markov chain), a\ is the variance in the vari- 
able x and p XJ is the correlation coefficient between variables 
x and y. The current values of the harmonic amplitudes Ay_ and 
Bk are determined as a function of the current values of the pa- 
rameters M u M 2 , 71, T 2 , R U R 2 using Eqns. (7), (1 1)— (19) and 
Kepler's third law. The medians, variances and correlation co- 
efficients are fixed to their values as calculated using the joint 
posterior distribution resulting from the earlier analysis of the 
light curve data (as described in §2). 

To ignore, in effect, the statistical influence of DB or ELV, 
we remove the terms in Xt associated with A\ [sin(</>)] or B2 
[cos(20)], respectively. In each scenario, a Markov chain of 
length 5 x 10 6 was generated and the resulting posterior dis- 
tributions were inspected to have adequately converged. 

We report the results of our analysis of the system parame- 
ters in Table [2] The two columns in the table are the system 
parameters determined from the ELV and ILL effects, on the 
one hand, and the DB and ILL effects on the other. We list 
the median values of the parameters, as well as the 15.4% and 
84.2% (corresponding to the ±l-sigma values) values of the 
cumulative distribution function for each parameter of inter- 
est, having marginalized over the remaining parameters. In 
the lower-right corner of Fig.|4] we show the joint posterior 
distributions in a number of parameters to show correlations 
and differences between emphasizing DB or ELV. 

Of the approximately dozen system parameters reported in 
Table 2 (luminosities, radii, T e ff, masses, gravity, and so forth), 
the agreement between the two columns (ELV and DB dom- 
inated analysis) is remarkably good except for those param- 
eters that involve the white dwarf mass. The latter depends 
heavily on the mass ratio which, in turn, is strongly depen- 
dent on the assumption of whether the ELV or DB amplitudes 
correctly reflect their respective physical quantities. 

The estimated masses, temperatures and radii of the sec- 
ondary, in either of the two scenarios, are consistent with a 
hot white dwarf companion. In §5, we discuss possible evolu- 



TABLE 2 

KHWD3 Model-Independent System Parameters 



Parameter 


ELV+ILL Model 


DB+ILL Model 


a [AU] 


06 1 -4-0 004 


061 -1-0 004 


1VI 2 I >Y1 1 


1 SO-l-0 006 


1 OS-l-0 00s 
u. lujznu.uuo 


Tr. /T 

I 2/ I l 


1 . JOIEU.U1 


1 .J JxU.Ul 


Mi [Mq\ 


Z.J ± U.J 


Z.J ± U.J 


f M L ,v 0J 


1.87±0.13 


1.86±0.12 


7j [K] 


9600±150 


9500±150 


L { [L©] 


26.2± 3.8 


24.8± 3.3 


Pi [g/cc] 


0.54±0.02 


0.55±0.02 


log(g!) [cgs] 


4.29±0.03 


4.30±0.03 


M 2 [Mq] 


0.37±0.08 


0.26±0.04 


Ri [«©] 


0.151±0.011 


0.150±0.010 


T 2 [K] 


14900±300 


14600±300 


L 2 [L ] 


1.00±0.15 


0.93±0.13 


Pi [g/cc] 


153±11 


109±9 


log(g 2 ) [cgs] 


5.65±0.04 


5.50±0.02 



Observed-Calculated Harmonic Amplitudes [x 10 





A, [sin(</>)] 


-0.48± 0.12 b 


O.OOi 0.06 




B[ [cos^)] 


0.01± 0.06 


O.OOi 0.06 




A 2 [sin(20)] 


-0.44± 0.04 b 


-0.44± 0.04 b 


(25) 


B 2 [cos(2<^)] 


-O.OOi 0.06 


-1.80± 0.35 b 


A 3 [sin(3<«] 


0.02± 0.04 


0.02± 0.04 




B 3 [cos(3</.)] 


0.08± 0.05 


0.02± 0.05 



NOTE. — a = Mass of the primary is assumed, b : 
the likelihood. See § 3 for details. 



Values did not affect 



tionary histories leading to these two outcomes and estimate 
the properties the system progenitors. 

3.6. Possible missing physics 

In the previous sections, we inferred system parameters for 
KHWD3 only after discounting certain aspects of our model 
for the out-of-eclipse harmonic content in order to resolve 
some discrepancies in the measured amplitudes that exceed 
the estimated statistical uncertainties. At this time we admit- 
tedly have no satisfactory explanation for these discrepancies, 
though it is comforting that they do not affect any of our sys- 
tem parameter determinations other than those involving the 
white-dwarf mass - and even then the white dwarf mass is 
uncertain by at most a factor of 2. We note that, to the best 
of our knowledge, this is the first time that all three effects, 
DB, ILL, and ELV, have been studied at such small ampli- 
tudes (e.g., parts in 10 4 ). However, in this regard, we would 
have surmised that the simple, classical, analytic models for 
these effects should work even better at these smaller ampli- 
tudes than at larger amplitudes. It is also possible that there 
are some additional, heretofore unexplained, physical effects 
which complement the harmonic content that we have not at- 
tempted to model, e.g. the O'Connell effect (see, for example, 
Davidge & Milone 1984). 

3.6.1. Non-zero eccentricity 

We measured a tentative non-zero value for the eccentric- 
ity e, at weak significance, with ecosOJ = 0.0029 ± 0.0005; 
however, we assumed a circular orbit in assessing the out-of- 
eclipse harmonic content. A small but non-zero eccentricity 
would lead to slightly modified expectations from all three ef- 
fects leading to the out-of-eclipse variation (ELV, DB & ILL). 
For very small eccentricity, its influence can be determined 
analytically with epicyclic approximations to the binary orbit, 
whereby the semi-major axis a/R\ in the model expressions 
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FIG. 4. — Joint posterior distributions (showing 68% and 95% confidence contours) and histograms for KHWD3 light curve model and system parameters. The 
upper left portion shows some results from the analysis, as described in §2, of the Kepler light curve data for KHWD3. See Table 1 for complete results. The 
lower right portion shows some results from a subsequent analysis, as described in §3, to determine absolute system parameters for KHWD3 given the results 
from the photometric analysis. The shaded contours and histograms correspond to the scenario in which only ellipsoidal light variation and illumination were 
assessed. The non-shaded contours correspond to the scenario in which only the Doppler boosting and illumination were assessed. See §3 for more details. 



(Eqns. 11-19) is replaced with the instantaneous separation 
r/R\, which is a function of <fr, ecostt! and esinES. We find 
that DB is the only effect that can have a significantly altered 
harmonic spectrum [to linear order in eccentricity and given 
our restriction on e cos 03]. This alteration is entirely mani- 
fested in the inclusion of a sin(20) mode whose amplitude 
relative to the dominant sin <fi mode is 



— wesinGJ. (26) 
A\ 

We did detect a non-zero amplitude A2/A1 = -0.43 ±0.06, 
however, this amplitude is much too large to be accounted for 
by eccentricity alone. 



3.6.2. Departure from synchrony 
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We have assumed synchronous rotation between the pri- 
mary, secondary, and orbit when calculating the contribution 
to the harmonic content from ELV or ILL (as described in 
§3.1). Should the primary be rotating at a different rate from 
the orbital frequency, frictional drag within the primary may 
force the tidal distortion due to the secondary to lag or lead 
the orbit (as is the case with the tides imposed on Earth by the 
Moon). The illumination effect may also lag or lead as a re- 
sult of both the modified phase variation of the sky -projected 
cross-sectional area and additionally, depending on the depth 
in the stellar atmosphere at which radiant energy from the sec- 
ondary is deposited, from a "hot spot" that is shifted from the 
substellar point. 

To investigate this possibility further, we carried out an 
analysis where we assumed that both the ELV and ILL ef- 
fects are shifted by a common variable phase from their nom- 
inal circulation. DB is independent of the assumption of 
synchrony and was not permitted a phase shift. Given this 
additional degree of freedom, we find that we may account 
for the measured values of all harmonic amplitudes includ- 
ing the relatively large, and previously unexplained, sin(20) 
mode. However, the results of this analysis remain unphysi- 
cal (Mi ss 0.8 Mq and R\ « 1.3 Rq). Moreover, we deter- 
mined that these effects lag the orbit by 5° while reasonable 
expectations for the efficiency of tidal dissipation (quantified 
via the parameter Q) would predict maximal phase offsets of 

1/(2 ~ 10" 6 . 

3.6.3. Other effects 

Here, we mention two other effects not included in our 
model that could significantly affect our interpretation of the 
harmonic content and/or the photometric results: 

• Stellar spotting - A persistent stellar spot on a syn- 
chronously rotating primary could induce variations 
with power in multiple harmonics, most likely in the 
fundamental mode (competing with DB & ILL). The 
amplitude of this effect depends on the size and flux 
contrast, relative to the unspotted photosphere, of the 
spot. While spots in tidally locked binaries may persist 
for months, it is unlikely that an A star lacking a sub- 
stantial convective outer envelope would exhibit signif- 
icant spotting behavior (Strassmeier 2009). A possible 
exception would be a magnetically active A star (e.g., a 
"peculiar" or Ap star; see Kochukhov (2010) and refer- 
ences therein). 

• Rapid rotation of the primary - In the previous section, 
we discussed possible phase lags in the ELV and ILL 
effects as a result of non-synchronous rotation. In ad- 
dition, the primary may be significantly oblate as a re- 
sult of very fast rotation. In this case, a/R\ as inferred 
in the light curve analysis may be different from the 
true normalized semi-major axis. In particular, if the 
spin axes of the rapidly rotating primary and the orbit 
are aligned, then a/R\ would likely be overestimated 
(Barnes 2009). Additionally, very fast rotation may in- 
duce an unexpected stellar brightness profile as a re- 
sult of strong gravity darkening; both the inferred ra- 
dius ratio, Rt/Ri, and a/R\ may be inaccurate as result 
(Barnes 2009). However, given the short orbital period 
of 3.3 days, it is most likely that the orbit and the pri- 
mary star have already synchronized (see, e.g., Torres, 
Andersen, & Gimenez 2010). 



3.7. Erroneous a/ R\? 

It should be noted that a smaller value for a/R\, closer to 
6.2 than its current estimate of 7 ±0.1, would resolve the dis- 
crepancy between ELV and DB amplitudes for Mi ~ 2.5 M Q . 
If this lower value where accurate, then q ~ 0.11, as would 
be estimated using the DB amplitude alone. However, the 
temperatures of the primary and secondary would be reduced, 
according to the ILL model, with lower values of a/R\. In 
particular we find, T\ ~ 8400 K for a/R\ = 6.15. Additionally, 
the inverse correlation between R 2 /R\ and a/R\ (see Fig. 4) 
suggests a relatively larger white dwarf radius for a smaller 
ajR x . 

We have taken care to correctly model the transit and 
eclipse light curves (see §2) but note that, given the extreme 
correlation between inclination and a/R\ (see Fig. 4), small 
errors, model inaccuracies, or non-uniform data sampling 
may lead to systematically inaccurate values for a/R\ that 
vary greatly in an absolute sense. 

We remark that the ^40 data points occurring during the 
ingress or egress phases of the secondary eclipse carry the 
overwhelming majority of the statistical weight in determin- 
ing a/Ri, i and ecos(ttJ). While the MCMC algorithm (§2.2) 
will likely yield a robust measure of parameter precision and 
correlation, the accuracy of the most-likely value is affected 
by correlated noise and the (sparse) distribution of data points 
during eclipse ingress or egress. Compiling future Kepler data 
for KHWD3 will most likely resolve any current statistical bi- 
ases. 

4. INFERRED MASS OF THE WHITE DWARF 

Recall from §3 that, depending on whether our analy- 
sis emphasized the ELV and illumination amplitudes or the 
Doppler boosting and illumination amplitudes, we find that 
q = 0.150 ± 0.006 or q = 0.105 ± 0.008, respectively (see §3 
for details). For each of these values of q, we would infer a 
white dwarf mass of 

M 2 = (0.37 ± 0.08) [ ^' ) M ELV & ILL (27) 
\2.5M / 

M 2 = (0.26 ± 0.04) ( — ! )Mq DB&ILL (28) 

(see Table 2). We discuss the consequences of these values for 
the binary stellar evolution of the system in the next section. 

5. SYSTEM PROGENITORS 

The presence of a white dwarf in KHWD3, coupled with 
the short orbital period (3.273 days), indicate that there was 
necessarily a phase of mass transfer/loss during the prior evo- 
lution of this binary. One obvious constraint on the primordial 
binary is that it must have had a total mass, M p +M s , equal 
to at least the current mass of the binary Mi +M 2 — 2.6M , 
where M p and M s are the initial masses of the primordial pri- 
mary (the white dwarf progenitor) and the secondary, respec- 
tively. As is readily evident, and we explore below, this im- 
plies that the mass of the primordial primary was greater than 
1.5M , and therefore had a radiative envelope at the time 
when mass transfer to the primordial secondary commenced. 
In turn, this indicates that the mass transfer from the primary 
to the secondary took place, at least initially, on a thermal 
timescale. 

Given the two different possible masses for the white dwarf 
(see §4), there are two slightly different evolutionary scenar- 
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ios for the formation of the white dwarf. We discuss each of 
these in turn. 

5.1. Primordial Primaries with M < 2.2M Q 

Even if the donor star is more massive than the accretor, the 
mass transfer can still proceed quite stably if the donor has 
a radiative envelope. The mass ratio of the primordial pro- 
genitor binary, q pmg = M p /M s , up to which the mass transfer 
is stable, depends on the masses and the orbital period when 
mass transfer commences, but can exceed q plog ~ 2. For stars 
with initial mass < 2.2M , a degenerate He core develops 
and there is a nearly unique relation between the radius of the 
evolving star and the core mass. This leads to a tight relation 
between the mass of the remnant white dwarf, once the en- 
velope of the primary has been completely transferred to the 
secondary and/or lost from the system, and the orbital period: 



P orb ~4.6x 10 6 



(l + 25M" + 29Mt ri ) 3 / 2 



days (29) 



where M w( j is in units of Mq (see, e.g., Rappaport et al. 1995; 
Ergma 1996; Tauris & Savonije 1999; Lin et al.2010). If 
we solve this non-linear equation for the value of M w d that 
matches the current orbital period of 3.273 days, we find: 



M wd = M 2 ~ 0.21 ± 0.02 M G 



(30) 



where the estimate in M w d is uncertain by — 10% (Rappaport 
et al. 1995). This mass range is marginally consistent with 
the lower of the two possible masses measured for the white 
dwarf in KHWD3 (0.26 ± 0.04M Q ). 

Another major issue that needs to be addressed regarding 
the white dwarf is how to explain its large radius (R2 = 0. 15 ± 
0.01 R Q ) and high effective temperature (T 2 = 14, 500 ±300 
K). In particular, the radius is approximately 7-9 times the 
radius of a degenerate He star of the same mass. According 
to cooling models of He white dwarfs formed in an analo- 
gous manner (i.e., after losing its envelope to a companion 
star; Hansen & Phinney 1998; Nelson, Dubeau, & MacCan- 
nell 2004; Ph. Podsiadlowski, private communication 2010) 
the radius can remain this bloated for a considerable period 
of time for very low-mass white dwarfs. The degree to which 
this is possible depends on the amount of residual H-rich at- 
mosphere the core retains after losing its envelope. In particu- 
lar, a He white dwarf of 0.21 Mq with a H-rich atmosphere of 
0.01 Mq can remain this large and hot for — 150 Myr, while a 
0.19 Mq He dwarf with a similar H-rich atmosphere can re- 
main thermally bloated and hot for up to 300 Myr. Even in 
the case of 0.21 Mq white dwarf, 150 Myr is a non-negligible 
fraction of the nuclear evolution time of the — 2.3M Q par- 
ent star. Degenerate He white dwarfs with masses > 0.25 Mq 
would not remain hot and thermally bloated for nearly long 
enough to have a plausible probability of catching them in 
such a state. Therefore, the lower limit of the measured range 
of 0.26 ± 0.04M Q would be strongly preferred in this sce- 
nario. 

5.2. Primordial Primaries with M > 2.2Mq 

Stars with mass > 2.2M© do not evolve with degenerate 
cores and their radius does not follow the radius-core mass 
relation discussed above. In a close binary system, such stars 
are less likely to produce a remnant He core with a mass 
as low as — 0. 19-0.22 Mq. However, they could produce a 
white dwarf that matches the more massive of our allowed 
solutions, namely — 0.37 ±0.08 M Q . Such non-degenerate He 



stars, once the envelope of the primordial primary has been re- 
moved, would be on the He-burning main sequence (as long 
as M core > 0.3 IMq), and they could be quite luminous. 

An analytic fitting function for the luminosity of naked He 
burning stars is: 



1.53 x 10 4 M, 



10.25 



-He : 



M 9 e + 29 .5M^ + 3 1 .2M h ile + 0.047 



(31) 



(Hurley, Pols, & Tout 2000; hereafter "HPT"). The observed 
luminosity of the white dwarf in KHWD3 is L 2 ~ 0.80 ± 
0. 13L Q . The He star mass that corresponds to this luminosity, 
according to the above fitting formula, is down near the end of 
the He-burning main sequence with Mn e < 0.32M Q . The ra- 
dius and effective temperature of such a He burning star would 
be -0.055# Q and -26,000 K, respectively (HPT). The cor- 
responding measured radius and effective temperature of the 
hot white dwarf in KHWD3 is about 2.5 times larger in radius, 
and about half as high in temperature. This could be nicely 
accommodated (and produce roughly the same luminosity) if 
the He star has a modest residual hydrogen-rich atmosphere 
(i.e., > O.O2M ; see Han et al. 2002; Han et al. 2003). The 
nuclear lifetime of such a He-burning star, i.e., near the end 
of the He-burning main sequence, is — 10 9 yr. Such a star 
would be at the low-mass end of what are known as subdwarf 
B (sdB) stars (see, e.g., Han et al. 2003, their figures 15 and 
17). 

5.3. Possibility of a Prior Common-Envelope Phase 

We note that unstable mass transfer leading to the current 
system probably could not have resulted in a common enve- 
lope ("CE") phase leading to the successful ejection of the 
CE. The ejection of a common envelope around such a low- 
mass core would have resulted in a post-CE orbital period that 
is considerably shorter than 3 days. The only way in which 
a 3-day post-CE binary would ensue is if the pre-CE orbital 
period were extremely long - in which case the white dwarf 
mass would be much higher than is observed. For an im- 
pressive study of such post common-envelope systems see the 
work by Parsons et al. (2010). 

5.4. Constraints on M p , M s and Porb.init 

We explore here what the range of possible and likely pa- 
rameters the primordial progenitor binary could have been. 
We know three independent binary parameters of the current 
system: M\, M2, and P or t,. The appropriate application of 
conservation of mass and angular momentum yield two con- 
straints - insufficient to uniquely identify, by themselves, the 
initial system parameters. However, there are two other con- 
straints, discussed here, that are sufficient to define a relatively 
narrow range in parameter space for the primordial binary. 

As mass is transferred from the primary to the secondary, 
a fraction of it, (3, will be retained by the secondary, and the 
remainder will be lost from the system. The ejected mass will 
carry away an average specific angular momentum which we 
denote as a which is in units of the specific angular momen- 
tum of the binary. Unfortunately, we do not know a priori 
either a or (3, but we can make an educated guess about the 
former. In terms of the quantities we have already defined 
earlier, we can write the mass retention fraction as: 



M x -M s 
M p -M 2 



(32) 
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Conservation of angular momentum then yields a relation be- 
tween the orbital periods before (P,) and after (Pf) the mass 
transfer phase: 



Pi=Pf 



M y +M 2 )\M 2 ) \Mi 



c 2 



where the powers, C\ and C 2 are defined as 

Ci=3a(l-/3)-3 
C 2 =-3a(l-j8)/j8-3 



(33) 



(34) 
(35) 



(see, Podsiadlowski et al. 1992; Rappaport, Podsiadlowski, & 
Horev 2009). In the scenario discussed above for forming a 
white dwarf of mass ~0.25 M , we can constrain P, to be 
shorter than Pf = 3.273. In the case where a more massive 
white dwarf formed from a primary star with M p > 2.2M©, 
the initial orbital period could, in principle have been longer 
than 3.273 days, and the orbit could have subsequently shrunk 
during mass transfer. However, for substantially longer initial 
P 01 b, the core mass would likely exceed ^0.4 M . 

Another constraint that we can impose is that the initial or- 
bital separation should be large enough so that not only can 
both the primordial primary and secondary fit within their re- 
spective Roche lobes, but that there is sufficient room for the 
primary to evolve a substantial He core before the mass trans- 
fer commences. This constraint can be written as: 



Rp(M P ) = = £/(?) [G(M P +M S )} 



1/3 



V 3 



(36) 



where £ is the fraction of its Roche lobe that is filled by the 
primordial primary of radius R p ; f(q) is the ratio of the Roche- 
lobe radius to orbital separation, a, which depends only on 
the mass ratio, q; and P, is the orbital period before the mass 
transfer commences. We take as a somewhat arbitrary but 
quite reasonable constraint: £ < 1/2. 

A final constraint is the ratio of the thermal and nuclear 
timescales of the primary to those of the secondary. We have 
shown (Rappaport et al. 2009) that these both scale roughly as 



Mp 



-2.6 



(37) 



We require that this ratio be smaller than ^0.6 so that, while 
the original primary evolved, the secondary (the current pri- 
mary) did not evolve substantially off the main sequence by 
the current epoch. On the other hand, we simultaneously re- 
quire that this ratio be larger than about 0.2 so that the thermal 
timescales of the two primordial stars not be too dissimilar. If 
they were very different, then it is likely that the mass transfer 
would have become dynamically unstable. In summary, we 
require that 0.2 < t\ /t 2 < 0.6. 

We adopt a nominal value for the specific angular momen- 
tum loss parameter, a = 1.0; however, we also consider values 
of a = 0.75 and 1.5. This is a reasonable set of bounds since, 
for example, the specific angular momentum parameter, a, at 
the L2 point where matter is likely to escape from the binary 
system, ranges from 1 .08 to 1 .56 for a mass ratio covering the 
rather wide range of 0.1 < M s /M p < 2. 

When we combine the constraints imposed by conservation 
of mass, angular momentum, ratio of timescales, and initial 
Roche-lobe filling factors we find an acceptable region in ini- 
tial parameter space that is shown in Fig. [5] 




M p 

FIG. 5. — Regions in the plane of Mp and M s where the progenitor binary 
that produced KHWD3 could have originated. The three shaded regions cor- 
respond to the angular momentum loss parameters a = 0.75, 1.0, and 1.50, 
from top right to lower left. The upper left and lower right boundaries of each 
region are determined by tj/t2 = 0.6 and 0.2, respectively. The upper right 
boundaries are set by P 01 b = 3.273 days, while the lower left boundaries are 
set by allowing the primordial primary to fill no more than half of its Roche 
lobe. The contours are of constant mass retention fraction, /3. 
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Fig. 6. — Comparison of the four known systems with white dwarfs that 
are relics of stable Roche-lobe overflow. The plotted show the orbital pe- 
riod vs. the white dwarf mass for KHWD3, KOI-74, KOI-81, and Regu- 
lus. The heavy curve is the theoretically derived P m b(M w ^) relation, going 
as ~M^ d for the lowest mass white dwarfs, while the pair of dashed curves 
marks the estimated theoretical uncertainties (see text for details and Rappa- 
port et al. 1995). 
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5.5. Comparison ofKHWD3 with KOI-74, KOI-81, and 
Regulus 

We end our discussion of the prior evolutionary history of 
KHWD3 by comparing this system with three other white 
dwarfs in close binaries that also likely formed during a 
phase of stable Roche-lobe mass transfer. These are KOI- 
74, KOI-81, and Regulus with white dwarfs in orbital peri- 
ods of 5.3 days, 23 days, and 40 days, respectively (Rowe et 
al. 1010; van Kerkwijket al. 2010; Gies et al. 2008; Rappaport 
et al. 2009). The measured white dwarf masses in these sys- 
tems are 0.22 ±0.03, -0.3, and 0.30 ±0.02 M Q . Including 
the lower-mass solution for KHWD3 (i.e., 0.26 ± 0.04 M ), 
we find a modestly significant correlation between masses and 
increasing P or b. This correlation is plotted in Fig.[6]where the 
solid curve is the theoretically expected relation for progenitor 
stars of mass < 2.2M Q (see eq. 29 above). The white dwarf 
masses are not very well determined, and there is only a small 
spread in all their masses by a factor of —1.5. Nonetheless, 
given the extremely steep theoretical dependence of P rb on 
the mass of the white dwarf, the known systems are roughly 
consistent with this relation. More such systems, with even 
better determined white dwarf masses, will be needed in or- 
der to supplement the kind of information currently being pro- 
vided by radio pulsar systems (Rappaport et al. 1995; Tauris 
1998; Thorsett & Chakrabarty 1999). 

6. SUMMARY AND CONCLUSIONS 

In this work we have utilized public Kepler data to identify 
and study in detail a third hot white dwarf in a close binary 
orbit. The system is formally known as KIC 10657664, and 
we have given it the shorthand name "KHWD3". The system 
is comprised of a white dwarf with T e fs — 14,500 K orbiting 
an A star of mass —2.3 M Q . The white dwarf is extremely 
thermally bloated with a radius of O.15/?0, some 7-10 times 
larger than its degenerate radius. 

We have combined the transit and eclipse data as well as 
the low-amplitude, out-of-transit/eclipse light curve to infer 
all the system parameters. The large size of the white dwarf, 
coupled with its high effective temperature, produces rela- 
tively large eclipse and transit depths (—2% and 0.7%, respec- 
tively). The out-of-transit/eclipse periodic light curves are of 
much smaller amplitude (in the range of — 1 - 5 parts per 10 4 ) 
and are interpreted as being due to Doppler boosting, mutual 
illumination, and ellipsoidal light variations (due to tidal dis- 
tortions). The deduced system parameters are summarized in 
Table 2. 

The only ambiguity in the determination of the system pa- 
rameters lies in the mass of the white dwarf. Our analysis 
leads to two somewhat distinct, and currently unresolvable, 
possible solutions: 0.26 ± 0.04 M Q or 0.37 ± 0.08 M Q . The 
difficulty in distinguishing between these two possibilities for 
the white dwarf mass lies in the fact that both the Doppler ef- 
fect and the ELV both basically determine the mass ratio, and 
these two amplitudes in the light curve produce somewhat in- 
consistent results in this ratio (at the —2.8 a level). However, 
it is important to note that the remainder of the important sys- 
tem parameters are well determined, in spite of this particular 
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ambiguity. 

We have also briefly explored what the possible progenitors 
of this system might have been. The system almost certainly 
formed after a phase of mass transfer from the primordial pri- 
mary star (now the white dwarf) to the secondary star (now 
the A star primary of the system). If the primordial primary 
had a mass less than —2.2 M Q , its core mass would directly 
determine the final orbital period (3.273 days), which predicts 
the current white dwarf mass to be 0.21 ±0.02M Q , in modest 
agreement with the lower of our two solutions. Such a low- 
mass white dwarf (i.e., 0.19-0.22M Q ) could remain hot and 
bloated for a substantial fraction (i.e., > 10%) of the lifetime 
of the current primary A star. If the primordial primary had a 
mass > 2.2M©, then it is likely that our solution for a more 
massive white dwarf is the correct one. In this case, however, 
the white dwarf would almost certainly be undergoing nuclear 
burning at the current epoch. This, in turn, would imply that 
the mass would have to be low (— 0.31 -O.32M ), i.e., near 
the end of the He-burning main sequence, in order to explain 
a luminosity of < 1 L Q . 

Likely primordial stars of mass M p — 2.2M© and M s — 
1.4M Q (see Fig. [5} would imply a highly non-conservative 
phase of mass transfer with perhaps half the transferred mass 
being lost from the system. For this particular set of illus- 
trative primordial binary masses, the mass ratio is sufficiently 
large to possibly explain the high mass loss rate. We have 
compared KHWD3 with the two other hot white dwarfs dis- 
covered with Kepler as well as the Regulus system. The initial 
system masses for Regulus (see Rappaport et al. 2009) were 
M p — 2.1Mq and M s — 1.7M Q , which are seemingly not too 
different than for KHWD3. However, the final orbital period 
of the white dwarf in the Regulus system is 40 days, which is 
considerably longer than the 3.3-day period of KHWD3. The 
difference likely lies in the fact that in the Regulus system 
the mass transfer was considerably more conservative (lead- 
ing to a much wider orbit; see Eqn. 33). This is quite plausible 
given that the initial mass ratio of — 1 .2 for Regulus was sig- 
nificantly closer to unity than for KHWD3 which may have 
been closer to — 1 .6. 
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